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1. Introduction 

Since the discovery of graphene -a single layer of carbon atoms arranged in a honeycomb 
lattice - it was clear that this t ruly is a unique material system with an unprecede nted combi- 
nation of ph ysical properties iGeim & Novoselovl i2007n:lNovoselov et alJ i2004l): Novoselov 
& Jiang et al. j2005l) ; lNovoselov & GeimI j2005l) ; IZhang et alJ j2005l) ). Graph ene is the thinnest 
memb rane present in na ture -just one ato m thick- it is the strongest material jCastro Neto et alJ 
^2009^ ). it is tra nsparent iNair et al. ( 2008)) a nd it is a very eood conductor ( Novoselov & Jiang 
et al. t2005.) ; , Novoselov & Geim. (.2005.) ) with room temperature charge mobilities larger than 
the typical mobilities found in silicon. The significance played by this new material system is 
even more apparent when considering that graphene is the thinnest member of a larger family: 
the few-layer graphene materials (FLGs). Even though several physical properties are shared 
between graphene and its few-layers, recent theoretical and experimental advances demon- 
strate that each specific th i ckness of few-la y er graphene is a material with u nique physical 
properties jBaeetal.1 <2010l) ; lKim et al.l b009l) ; lLee etaD b008l) ; lNair et alj j2008l) ). 
All few layers graphene are good conductors. However, striking differences in the nature of 
this conductive state emerge when a perpendicular electric field generated by gate voltages 
is applied onto the few-layers. In a single layer graphene transistor, the current is modulated 
by a gate voltage but it cannot be switched off since in the energy dispersion of graphene 
ther e is no band-eap (valence and conduction bands touch each other) ( Novoselov & Jiang 
et al. i2005l) ; lNovoselov & GeimI i200j)). Recent experimental advances showed that bilayer 
graphene, also characterized by touching valence and conducti on bands, develops an energy 
gap when subjected to an external perpendicular electric field iZhang et alJ ^20091) ). Bilayer 
graphene is the only known semiconductor with a gate tuneable band-gap. Opposed to the 
case of single- and bi-layer, the trilayer material is a semimetal with a gate-tuneable band 
overlap between the conduction and the valence band. Indeed, the co nductivity of trilaye rs 
increases when a perpendicular electric field is applied onto the system jCraciun et al.l i2009l) ). 
The variety of physical properties found in different FLGs is the true strength of these newly 



discovered materials, and yet very little is known on few-layer graphene with more than 3 
layers. 

2. Electronic properties of few-layer graphene materials 

The electronic properties of a material are intimately related to its energy dispersion. There 
are several approaches to calculate the electronic energy bands and here we review the cur- 
rent understanding of th e graphene ma terials band structure within the non-interacting tight- 
binding approximation iWallacd Jl947ll ). When two or more carbon atoms are brought to- 
gether to form a regular lattice -such as the hexagonal lattice for a single layer graphene- the 
valence electrons of the different atoms interact. This leads to a broadening of the electronic 
eigenstates and ultimately to the formation of the continuous bands of a solid. 
An isolated carbon atom has 6 electrons ls-^2s^2p^, where the energies of the s-orbital and p- 
orbitals of the second electronic shell are very similar. Consequently, carbon can form a num- 
ber of hybridized atomic orbitals characterizing different geometries. In the case of graphene, 
one s-orbital and two p-orbitals {px and py) undergo a sp2 hybridization with a characteristic 
planar trigonal symmetry with an angle of 120°between each bond. This is the reason why 
each carbon atom within graphene has three nearest neighbors at a distance of = 0.142nm. 
On the other hand, the pz-orbitals overlap sideways with regions of highest electron density 
above and below the graphene plane and the energy dispersion of these /r-bonds determines 
the electronic transport properties of graphene materials. 




Fig. 1. Panel (a) shows the crystal structure of monolayer graphene whose imit cell contains 
two equivalent carbon atoms -A and B. The 3D plot in (b) shows the energy dispersion of 
graphene (see Eq. The valence and conduction band touch in 6 points, known as valleys. 
Whenever the onsite energy symmetry between the A and B sublattices is broken a band-gap 
opens in the energy dispersion of graphene as shown in Panel (d) -energy dispersion obtained 
considering H/^y^ ^ Hbb with H/^y^ = O.I70 and H^b = — 0.l7o- 

The hexagonal lattice of graphene is a composite lattice with two carbon atoms in the unit cell 
-indicated by A and B, see Fig. [1^- and basis vectors: 



ai = -V3/2aQi+3/2ao] 
a2 = V3/2aot +3/2ao]. 



(1) 



We consider the linear combination of atomic orbitals (LCAO) of the Bloch wavefunctions 
corresponding to the sublattice carbon atoms A and B of the form: 



= ;^E„ <P,(r-d,«), i = A,B. (2) 

Where N is the number of cells considered, k is the wave vector, d,„ is a lattice translation 
vector and (p, (r — d/„) is the pz local atomic orbital. The /r-orbitals electronic band structure is 
therefore a solution of the Schrodinger equation: 

H<f(k) = E(k),f (k) (3) 

where U is the Hamiltonian for an electron in the atomic potential given by the atoms in the 
graphene lattice and i^(k) is a linear combination of Bloch wavefimctions which for simplicity 
we can write in the form (f = a(^j^^ + b</)B, where (\>a and are given by Eq. |2]and a and h 
are two coefficients. To solve Eq. |3]we need to find the matrix elements Ef,y =< (pi\H\(pj > 
of the Hamiltonian and the overlaps between the Bloch wavefunctions < (pilcpj >■ We start 
from noticing that if the two carbon atoms forming the graphene sublattice are energetically 
equivalent, the onsite energies of the sublattice A and B are equivalent (Haa = Hgg) and 
without loss of generality we can set this energy equal to zero. The solution of Eq|3]is further 
simplified if we consider that the mo st significant hopping parameter is the first neighbour 
(7o « 2.8eV, iCastro Neto et al.l HqO^)) and that HgA is simply the complex conjugate of Ef^g 
(Hba = Hab*)- Therefore we can calculate the integral Ef^g =< (/)/i|Ef|(fg > where each atom 
is surrounded by three neighbours with relative coordinates (1/3, 1/3), (1/3,-2/3), (-2/3, 1/3). 
The term Ef^g reads: 

We can now project the solution of Eq. |3]onto < (pj\ \ and < if g | to obtain the system: 
aH^A + bHAB = E(k)a / foEf^g = E(k)fl 



uHba + bHBB = E(k)fo ^ \ aHAB* = E(k)b 

which has non-zero solutions for the coefficients a and b only if its secular determinant is zero. 
This condition leads to the energy dispersion of the graphene zr-orbitals: 



E(fc,.,/cy) = ±7oyi+4cos(^fcy)cos(^fc0+4cos2(^fc,) (6) 

where r = floV'S- The energy distribution vanishes at six points in the reciprocal lattice space 
with coordinates ±27r/r(l/v^, 1/3), ±27r/r(0,2/3), ±27r/r(-l/v^, 1/3), see Fig. [ijj. In 
these six K-space points the valence and conduction band touch one another, but only two of 
these points are independent. These are commonly indicated by K and K' and also known as 
valleys. The electronic states close to the Fermi level (E=0) are readily described by a Taylor 
expansion of the energy dispersion in Eq. |6] at a chosen K point. This reveals that charge 
carriers in graphene are mass-less Dirac electrons obeying a linear energy dispersion: 

E(fc) = ±i;F|fc| (7) 

with vr = \/3'Ynr / 2h the energy independent Fermi yelocitv jliang et alj i2007l): Novoselov & 
Geim j2005|) : IZhangetalJ<200g) ; IZhang&BraretalJ<2008|) ). 



The absence of a band-gap in the energy dispersion of graphene implies that the conduction 
in this material cannot be simply switched on or off by means of a gate voltage which acts 
on the position of the Fermi level, limiting the use of graphene in conventional transistor 
applications. Indeed, even when the Fermi level in graphene devices is at £ = 0, the current in 
graphene is far from beeing completely pinched-off. However, the gapless energy dispersion 
of graphene is a consequence of the assumption that the electron onsite energy between the 
A and B sublattice carbon atoms are equal. Whenever Hyi^ ^ Hgg a band-gap opens in the 
energy spectrum of graphene, see Fig. [l]:. A viable way to experimentally engineer such a 
band-gap consists in growing and /or depositing graphene on a commensurate honeycomb 
lattice formed by chemically inequivalent atoms which ultimately will originate a difference 
in the onsite energy between the A and B sublattices. So far, hexagonal BN is considered 
to be one of the most promissing candidates for graphene band-gap engineering since it has 
an almost commensurate c rystal structure t o the one of graphene and it has two different 
elements in each sublattice jGivannetti et al.l bOOTi) '). However, despite the growing interest in 
graphene on h-BN, no experirnenta l evidence has been reported yet of a band-gap opening in 
graphene on h-BN jXue et al.1 j201ll) ). 

On the other hand, bilayer graphene offers a unique alternative to the problem of band-gap 
engineering for reasons w hich will become clear when considering the bilayer non-interacting 
tight binding description jMcCanr] j2006l) ). The most common stacking of graphite planes 
found in nature is of Bernal type, where the A atoms in one layer are alligned on top of B 
atoms of an adjacent layer. The unit cell of bilayer graphene consists of a basis of four atoms 
labelled A\, Bl, AT. and B2 belonging to different atomic planes as indicated by the numerical 
index, see Fig. and b. 

Similarly to the case of a single layer we adopt the LCAO method with Bloch wavefunctions 
corresponding to the sublattice carbon atoms of the form: = e^'^<« cp j{i — d,„) with 

i = Al, Bl, A2, B2 and d,„ the translation vector of the i sublattice. For simplicity, we start 
by considering non-zero only the nearest neighbour coupling -i.e. 70 hopping from Al to Bl 
and fro m A2 to B2 atomic sites- and 7^ « 0.39eV the interlayer coupling between A2 and Bl 
atoms jCastro Neto et aP (|2009|) '). We also assume that all the carbon atoms lattice sites are 
energetically equivalent -i.e. Hy\jAj = H^jBj with i,] = 1, 2. In this case we consider the linear 
combination of Bloch wavefunctions of the form (p = acpy^i + bcpgi + c(p/^2 + d<pB2, with a, h, c 
and d coefficients. The solution of the Schrodinger equation is readily obtained projecting the 
solutions of Eq. |3]onto the states < (pi \ with i = Al, Bl, A2, and B2. The equivalent system to 
Eq. |5]for the case of bilayer contains 4 equations, and imposing the condition that the secular 
determinant is zero leads to the set of the bilayer four bands: 



E{k) = ±f ± ^ f + (hvphf (8) 

with i^F = \/37oJ'/2^. The energy dispersion in Eq. |8]showes that the lowest energy con- 
duction and valence band of bilayer graphene touch each other in each of the two K and K' 
valleys. The higher energy bands are instead shifted ±0.4 eV away from the £ = Fermi level 
position. In the low energy limit, the energy dispersion in bilayers is parabolic (see Fig. [2]:): 



E{k) = ±hV/2m*, 
where m* is the charge particles effective mass m* = 7i/2i;|. 



(9) 



Fig. 2. Panel (a) and (b) show respectively a top view and a 3D view of the crystal struc- 
ture of AB-stacked bilayer graphene whose unit cell contains four equivalent carbon atoms 
-Ai, Bi,A2 and B2. The hopping parameters 70, 7i and 73 are highlighted in the scheme. The 
graph in (c) shows the energy dispersion of bilayer graphene without an energetical asym- 
metry between the sublat tices (con t inuou s line) and with various values of the onsite energy 
difference (dashed lines) dKoshinol b009l) ). When all the hopping parameters are considered 
in the tight binding calculation, the low energy dispersion of bilayers is not any more simply 
parabolic but 4 Dirac cones develop. Panel (d) shows an equi-energy contour plots of the low- 
est electron band of bilayer graphen e for equi-energy onsite simmetry case with the 4 Dirac 
cones clearly visible jKoshinol ilOO^ ). 



If we lift the energy degeneracy between the carbon sublattices on t he two layers the gap-les s 
energy dispersion of bilayers develops an energy gap around E = (Castro et al. (2007; '201 Ol): 
[ kuzm e nko et al. ( 2009); Mak et al. (200 9); McCarm (2006); McCann et al. (2006); Oht a et M 
T2006h:IOostinea et a l. (2008); iRusso et all (l2009l);lTavc hatanapa t et all (l2010l);IXia et all ( I2016I) ; 
IZhang & Li et alj (|2008t) : .Zhang et all (|2009t) ; IZhou et al., C2008,) : ,Zou et aljj2010|) l The origin of 
this band-gap is readily imderstood when considering an onsite energy asymmetry A between 
the sublattices Al and B2. In this case the energy dispersion opens an energy gap equal to 2A 
at each valley in the energy dispersion. This condition is easy to realize experimentally in 
double gated de vices (see Section [3) and it led to the discovery of a tuneable band-gap in 
bilayer graphene dCraciun et al.l (|2010|) 'l. 

Recent advances in ultraclean suspended graphene devices made it possible to study details of 
the energy dispersion of bilayer in the vicinity of E=0 and at energy scales much smaller than 
the nearest neighbour hopping, see Fig. |3] In this case, the tight-binding description which 
includes higher order hopping parameters reveals that the energy dispersion of bil ayers has 
a trigonally deformed dispe r sion (see Fig. 3i) wi th 4 touching points in each valley teoshinj 
j200^ : lMcCann et alJ <2006|) : lMikitik et al .1720081) ). In the vicinity of each touching point the 
energy dispersion is linear (in the low energy 4 Dirac cones appear) and not parabolic as 
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Fig. 3. Panel (a) shows a scanning electron microscope image at an angle of 60°of a suspended 
2-terminal bilayer graphene device, obtained by wet-etching of the Si02 substrate. The graph 
in (b) is a low temperature (T=0.3K) electrical transport measurement of the suspended bilayer 
(charge mobility of « lOOOOc m^ /Vs). A series r esistance of 1.2KOhm has been subtracted to 
account for contact resistance iRusso et alJ i2010|) ) and for the resistance of the electric lines of 
the cryostat. 



predicted within the nearest neighbour hopping approximation. Therefore, depending on the 
energy scale, the Dirac fermions of bilayer graphene can lose or gain an effective mass. 
The topological discontinuity of the Fermi surface which occ urs when crossing from the Dirac 
cones to parabolic energy dispersion is a Lifshitz transition iLifshitzl il960|) ). It is rather ob- 
vious that a discontinuous change in the topology of the Fermi surface will be reflected in 
singularities of the thermodynamic and kinetic observables of the system. Since in bilayer 
graphene transistor devices the Fermi level can be continuously tuned by a gate voltage, this 
system is ideal t o study the occurrenc e of the Lishitz transition and its implication on the trans- 
port properties jLemonik et al.l j2010l) ). Currently, fundamental questions such as the stability 
of these 4 Dirac cones against electron-electron interactions and/ or mechanical deformations 
are at the focus of both theoretical and experimental research. 

For multilayers of more than two layers, there are two known structures called ABA (hexag- 
onal or Bernal) and ABC (rhombohedral) with different stacking manners as shown in Fig. 

and b. For ABA multilayer graphene, the Hamiltonian can be approx imately decomposed 
into a superposition of the monolayer-like and bilayer-like subsystems jKoshino et alJi200^ 
I2OO8)). Specifically, 2M-layer graphene contains M bilayer-like subbands, while (2M -|- 1)- 
layer graphene has M bilayer-like subbands and one monolayer-like sub bands. The subband s 
appearing in N-layer grap hene are characterized by a single parameter ^Guinea et al.l i2006l) ; 
iKoshino" et al.1 l|2007l;l200^ ): 



^ = 2^°^ NTT' ^''^ 

where fj = 1,2,--- ,[(N-|-l)/2] is the subband index and [(N-|-l)/2] is the maximum inte- 
ger which does not exceed (N + 1) / 2. The sub-Hamiltonian of A 7^ is equiavelent to that of 
bilayer graphene with the nearest interlayer coupling parameters multiplied by A, while that 




a) ^) 
Fig. 4. Panels (a) and (b) show respectively the crystal structure of Bernal stacked and rhom- 
bohedral stacked trilayer graphene with the complete set of tight binding parameters. 



of A = 0, appearing in odd N only, becomes that of monolayer graphene. For instance, the 
trilayer-graphene comprises six bands of which two have linear (monolayer-like) dispersion 
and four have parabolic disperion (bilayer-like). On the other hand, ABC multilayers have 
quite different electronic structures: the low-ener gy spectrum contains only a pair of conduc - 
tion and valence bands touchin g at zero energy l|Guinea et al ] j2006|) : lKoshino et all l|2009t) ; 
iKoshinol l2oToh : lMin et alj bOOSi) ). These two bands are given by the surface states localized at 
outer-most layers, and the dispersion is roughly proportional to , and becomes flatter as N 
increases. 
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Fig. 5. The graphs show the self-consistently calculated band structures in trilayer graphene 
near the K point respectively in the first neighbour approxini ation (i.e., 70 7i only) and in 
the full parameter model, adapted from lKoshino et all j2009l) . 



The response of ABA multilayer graphene to the gate electric field perpe ndicular to the layer s 
is unique and qua litatively different from that of AB-bilayer graphene dCraciun et alj ( |2009|) ; 
IRusso et alj j2009|) '). In trilayer graphene, for example, the gate field breaks mirror reflection 
symmetry with respect to the central layer and hybridizes the li near and p arabolic low-energy 
ban ds, leayine the semi-metallic band overlap near zero enerey jCraciun et al. (2009 ): Koshino 
et al. i200<^ ). Only the tight binding model based on the full set of hopping integrals (i.e. 70 
to see Fig. \5\ can succesfuUy describe the experimentally observed band-overlap ( Craciun 
et al. j2009l) :l Russo et alj j2009l) . The tuneable semimetallic nature of trilayers together with 
the tuneable band-gap in bilayers demonstrates that graphene materials are characterized by 
a unique range of physical properties not found in any other known material system. Theo- 
retically, a similar band overlapping is expected in thicker ABA multilayers as well, though 
very little is known experimentally. 

Finally, graphene multilayers can also be constructed by carefully transferring layer by layer 
individual graphene sheets on a substrate. The rise of this new family of materials, the meta- 
few layers (picTa = "beyond" the few layer materials), offers the unique possibility to control 
the hopping parameters by stacking engineering. In particular, the interlayer hopping param- 
eters are responsible for the rich low-energy band dispersion of these materials. Therefore, 
the ability to control these hopping parameters by means of graphene stacking engineering 
-for instance making thin tunnel barriers between two subsequent graphene layers- holds 
the promise for unprecedented fimctionalities as compared to natural few-layers. This layer- 
by-layer engineering of graphene materials has already been successfuly employed to create 
transparent electrodes in organic solar cells, and an efficiency similar to the ITO electrodes 
was recently demonstrated ( Wang et al. (201l|)). 

The physics of graphene has expanded at a rapid pace mainly thanks to the easily accessible 
electronic properties in simple transistor geometries, see Fig. |6^. In these devices, metallic 
contacts inject charge carriers into the conductive graphene channel, whereas the Fermi level 
of graphene is continuously driven from the valence to the conduction band by means of a gate 
voltage. As the Fermi level is driven inside the conduction (valence) band, the conductivity 
increases with increasing the concentration of electrons (holes) induced by the gate voltage. 
At the touching point between the valence and conduction bands the Fermi level crosses the 
zero density of states point -i.e. the Dirac point- where the conductance reaches its minimum 
value. Indeed, despite the density of sates vanishes at the Dirac point the conductivity re- 
mains fini te with a minimurn value of w 4e^ / nh for the b allistic t ransport regime (G eim & 
Novoselov l2007l) ; lKatsnelsonl j2006l) ; lNovoselov & GeimI j2005); Tw orzydlo et alj <2006l) ). The- 
oretically a finite conductivity is expected for ballistic electrical transport in graphene at the 
neutrality point assisted by evanescent wave propagation. This evanescent wave propagation 
gives rise to a unique non-monotonous dependence of the Fano factor on th e charge density 
in shot-noise jPanneau & Wu et al.l bOOSi) ; iDanneau & Craciun et al.l i2008l) ; iDanneau et al.l 
<2009|) : lTworzydlo et alJ (|2006h ). The predicted minimum conductivity 4e^ / nh has only been 
experimentally observed in devices with a short graphene channel. So far, most of the experi- 
mental studies on the minimum of conductivity have been conducted in supported graphene 
-e.g. graphene on Si02 substrate. Ultrapure suspended graphene devices do not suffer of 
the presence of substrate-related disorder and are therefore the ideal candidate to address the 
physics governing electrical transport at the minimum of conductivity. 

Indeed, the presence of disorder such as adatoms or molecules, ripples of the graphene sheet, 
atomic vacancies and topological defects are expected to affect the electronic properties of 
graphene such as the conductivity and charge carrier mobility. In particular, close to the Dirac 



point, charged impuriti es create elec t ron/h ole puddles which domir\ate the charge transport 
properties of graphene iMartin et alj i2008ll ). For small energies around the Dirac point recent 
theoretical advances also pointed out that strong short-range interaction caused by scattering 
off adatoms and molecules on the surface of graphene can actually be the dominant source of 
disorder limiting the charge carrier mobility. During the fabrication process of transistor de- 
vices, graphene is exposed to environmental conditions, it is therefore likely that for instance 
hydrocarbons are adsorbed on the surface of graphene. Whenever an hydrocarbon covalently 
bonds to the graphene, the 2j>z delocalized electrons are localized into a cr-bond, i.e. a co- 
valently bond hydrocarbon effectively act as a vacancy. The charge carrier scattering off the 
resonant states induced in the vicinity of the Dirac point together with scattering off charged 
impurities are likely to play a dominant role in gra phene devices, and this is current ly at the 
focus of both theoretical and experimental research iFerreira et al.l(|201l|) ; rPerejj2010|l l 

3. Experimental observation of gate tuneable band structure in few-layer graphene 

One of the most remarkable physical property of graphene materials is the ability to reversibly 
tune the band structure of these svstems simplv bv means of an external electric field ( Craciun 
et al. |2010|)). In standard semiconducting materials a precise value of the band-gap is engi- 
neered during the growth process, therefore the value of this energy gap carmot be reversibly 
controlled in situ in a device. Few-layer graphene is the only known class of materials to 
exhibit a gate tuneable band stru cture and this unpre cedented property paves the way for 
devices with novel functionalities (jCraciun et all l|2010|) ). 

Different experimental approaches have been implemented to address the ba nd structure 
of graphene materials. Respectively, pioneering charge transport experiments dCastro et al.l 
12007); Craciun et al. 12009); Oostinga et al. 1 200 ^)) followed by photoemission spectroscop y 
Ibhta et all |200d) ) and infrared sp ectroscopy jMak et all | |2009|) ; IZhang & Li et all bOOSi) ; 
IZhang et alj i2009ll ; lzhou et al] ( l2008h ) have highlighted complementary aspects of the energy 
dispersion of few-layer graphene with different number of layers. 

Possibly the best device geometry in which to address the electric field tuneability of the FLGs 
energy dispersion is a double gated design, where the graphene materials are sandwitched be- 
tween a top- and a back-gate, see Fig. |6^. Double gated geometries have a dual valency, that 
is they offer a simple way to independently and contin uously control in-situ both the band 
structure and the Fermi level by means of gate voltages dCraciun et al. I j2010ll ). A typical de- 
vice layout comprises a source and a drain contact to a graphene flake exfoliated onto Si02/Si 
substrate -which serves as a back gate- and a nanofabricated top gate. Since the conductance 
is determined by the features of the energy bands in the thermal shell kg T around the Fermi 
energy (ep), any modification of the conductivity in response to the perpendicular electric 
field is purely a consequence of the changes in the energy dispersion. This intimate relation 
between conductivity and energy dispersion has allowed the dis covery that b ilayer graphene 
is the only known semiconductor with a gate tuneable band-gap jZhang et al.l (2002)) ^rid con- 
versely that trilayer graphe ne is the only knovm semimetal with a gate tuneable conduction 
and valence bands overlap dCraciim et alj ( |2009|) ). 

A perpendicular electric field applied onto the few-layer graphene materials breaks the en- 
ergetic symmetry between the planes of these multilayer systems (see Section |2). This asym- 
metry is then reflected in the energetic inequivalence between carbon atoms belonging to dif- 
ferent sublattices -which in multilayer graphene belong to different layers. Experimentally it 
was observed that the in-plane electrical transport properties of each specific FLG thickness 
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Fig. 6. (a) Schematic picture of a double gated bilayer graphene device, (b) Square resistance 
Rsq of bilayer graphene as a function of back gate voltage (V^q) measured for different fixed 
values of the top gate voltage (Vjq) at T=300mK. The position of the Fermi le vel and the 
applie d perpendicular electric field are controlled by VgQ and Vjq, adapted from lRusso et al.l 
i2009|) . (c) the square resistance of a bilayer device versus Vf for different values of the back 
gate voltage Vi, and (d) shows the corresponding infrared abosorption spectra, adapted from 
|Zhang et al.' 1 2009). D represents the average value between the top- and back-gate electric 
displacements applied to the bilayer. 



change in a unique way in response to a finite external electric field (Ecx)- In all cases the re- 
sistance exhibits a maximum whose value and position in gate voltage depend on the voltage 
applied to the gate on which a fixed potential is applied during the measurement. 
In bilayers whenever E^x 7^ V m^^ the maximum of square resistance (R™^) displays a 
metal-to-insulator transition consequence of the opening of an electric field induced gap in 
the energy dispersion. Once the energetic equivalence between the sublattices is restored (i.e. 
Eex = 0), the energy gap reduces to zero. Consistently, charge transport experiments have 
reported large on/ off ratios of the current in double gated graphene bilayers when the Fermi 
level crosses from the conduction (valence) band through the band-gap, see Fig. ^ and c. 
However, the values estimated for the band-gap from transport experiments are systemati- 
cally much smaller than the theoretically predicted energy gap. Typically in transport a mo- 
bility gap is observed on an energy scale of a few meV for an average perpendicular electric 



displacement field of D = TV I nm, whereas a band gap of 200 meV is theoretically expected. 
Furthermore, the specific temperature dependence of J?™^ measured in bilayer graphene is 
incompatible with a simple thermally activated transport over a band gap, but it exhibits the 
functi onal dependen ce typic al of varia ble range hopping due to a finite sub-gap density of 
states jOostinga et al. (2008); iRusso et al. (2009); Taychatanapat et alj j2010l) ; lYan et a l. (2010i); 
IZou et al.H2o'lo|v i 

The dichotomy of a mobility gap in transport experiments and the theoretically expected en- 
ergy gap is currently fuelling both theoretical and experimental discussions. Though sev- 
eral possible explanations have been put forward, transport studies in different geometries 
demonstrate that the temperature dependence of the conductance in bilayer graphene in the 
diffusive regime can be explained by the parallel of a thermally activated current over the 
energy-gap and a variable range hopp i ng cu r rent through the disorder induced sub-gap states 
jTaychatanapat et al.l f20Kl) : lYan et atl ( l2010l) ; IZou et al.ll2O10l) ). 

The first direct observatio n of a gate-tuneable energy gap in bilayers was reported in in frared 
spec troscopy experiments llMak et al J ll2009h :IZhang & Li et al.1 ll200g) :IZhang et all <2009h •Zhou 
et al. boo g)). This technique is mostly sensitive to band to band transitions, therefore it is not 
affected significantly by transitions in the sub-gap energy range. Whenever a band-gap is 
open in bilayer graphene, the infrared absorption displays a highly intense peak in the ab- 
sorption spectra corresponding to the transition of charge particles from the top of the valence 
band to the bottom of the conduction band, see Fig. |6}l. This peak in the infrared absorption 
spectra corresponds to the energy gap between the valence and conduction band and it has a 
pronounced gate tuneability. In this way an electric field induced gap ranging from meV up 
to 250 meV at D = 3V / nm has been reported, which is consistent with theoretical predictions. 
At the same time, transport experiments in the very same devices show a small increase of the 
maximum resistance as a function of Eg^ confirming the presence of a large disorder induced 
sub-gap density of states. 




a) 




Fig. 7. (a) shows the Hall coefficient as a fimction of back gate voltage V^q (open circles) for a 
fixed perpendicular magnetic field of 9T at 50mK for a trilayer graphene device. The red curve 
is a fit. The insets depict schematically the position of the Fermi level (ef ) at different values 
of VgQ. The graph in (b) shows a 4-terminal resistance measurement of the trilayer device in 
the inset of (c) versus top gate and for different values of back gate, (c) shows the electric field 
dependence of the b and overlap Se for the same trilayer device. All panels are adapted from 
ICraciunetaIH2009l) 



Contrary to bilayers, ABA -stacke d trilayers d i splay a decrease of Rfq"^ with increasing E^x, 
see Fig. jCraciun et al.l j2009h : lRusso et aD ^0091) ). The overall electric field dependence 
of the resistance of trilayers can be explained adopting a two band model with an electric 



field tuneable band overlap between the conduction and valence band (Se). To this end tri- 
layer graphene is the only known semimetal with a gate tuneable band overlap, see[7]: . This 
unique property was independently demonstrated by magneto-transport experiments of the 
Hall coefficient (Rh), see Fig. [7^. In particular, measured at a fixed external perpendicular 
magnetic field displays a characteristic sign reversal corresponding to the cross-over between 
different types of charge carriers involved in the conduction (electrons and holes). The back 
gate voltage range over which Rf{ changes sign gives a band overlap w 28 meV, see Fig. [7^. 
To date little is known experimentally on thicker few layer graphene materials with more than 
3 layers. Recent infrared spectroscopy experiments address the evolution of the electronic 
properties from the one of mass-less Dirac electrons in a single layer to the massive particles 
of bulk graphite, presenting a systematic study from 1 up to 8 Bernal stacked graphene layers. 
Measurements of infrared conductivity show that the key features of the 2D band structure 
of few-laye r grap hene can be achieved on the basis of zone folding of the 3D graphite bands 
jMak et al.1 bOld) ). However, so far the electron transport properties of these thicker few- 
layer graphene materials is largely unexplored, preventing us from identifying the best suited 
thickness of few-layer graphene for a given application. 

4. Landau level structures in few-layer graphene 

When the charge carriers travelling in few-layer graphene experience a perpendicular mag- 
netic field (B), their trajectories are bent due to the Lorentz force. In the quantum regime, 
these cyclotron orbits give rise to discrete energy levels known as Landau levels. The precise 
sequence of these Landau levels reflects the nature of the charge carriers in the few-layer 
graphene under consideration. In particular, the Landau levels sequence for single layer 
graphene is intimately related to the mass-less nature of the graphene Dirac fermions (see 
Eq. and it is very different from what is known in conventional two-dimensional electron 
gases. 

The Landau level energies of monolayer graphenes are given by l |McClurel(|l956|) ): 



E„=hu)BSgn{n)J\n\ (w = 0, ±1, ±2, ..), (11) 



with fiLOg = sJlnv^eB. Each level is four fold degenerate, that is spin and valley degener- 
ate. Due to the linear energy dispersion (see Eq. |7), the energy spacing between the Landau 
levels is proportional to \/B rather than B as in usual two dimensional systems. At a fixed 
value of external magnetic field the energy gap between Landau levels in graphene are much 
larger than the corresponding gaps opened in other 2DEGs (for B=1T the energy gap between 
n=0 and n=l in graphene is approximately 35 meV). Another unique feature of the graphene 
Landau level spectrum is the existence of a zero energy level. This causes the half-integer 
quantization of Hall conductivity per spin and valley, and is also responsible for the huge 
diamagnetic susceptibility characteristic of this system. Specifically, the orbital susceptibility 
has a strong singularity at band touching point (Dirac point), which at zero temperature is 
expressed as a function of Fermi energy £f as: 

2 2 

X{^f) = -gvgs^Siep), (12) 



where gv = gs = 2. are respectively the spin and valley degeneracies <McClurell|l956h ). 



The low -energy Landau leve ls dispersion of bilayer graphene is approximately given by the 
relation jMcCann et alj ( |2006|) ): 



Esn=shcOc^n(n-l), (s = ±, w = 0, 1,2, ..) (13) 

with cvc = eB/m* the cyclotron frequency associated with the effective mass m* of bilayer 
graphene. The Landau levels energy spacing is now linear in B owing to the usual quadratic 
energy dispersion of bilayers, see Fig. |8ji. The two lowest levels of n = (per spin and valley) 
appear at zero energy. This amounts to 8-fold degeneracy in total and causes doubling of the 
Hall conductivity jump at zero electrori dens i ty, see Fig. [8b. The orbital susceptibility for small 
Fermi energy becomes teoshino et al] j2007l) ; ISa"fran et all jl984l) ): 

which has a logarithmic singularlity in contrast to the delta-function in monolayer graphene. 
The sigularity is weaker than in monolayer, since the Landau level spacing is narrower so 
that the total energy gain in magnetic field at Ef = becomes smaller. When increasing the 
magnetic field amplitude, the energy of the particular Landau level other than zero-energy 
levels crosses over from linear B to \/B , in accordance with the crossover of the zero-field 
dispersion from linear to quadratic. 



dG/dB 




a) B(T) b) Vg(V) 



Fig. 8. Panel (a) is a color coded plot of dG/dB versus gate voltage for the magnetic field 
range of OT to 2T and at T=0.3K measured in the suspended bilayer device of Fig. |3] As the 
magnetic field is increased in addition to an insulating state at filling factor v = 0, also the 
quantum Hall plateaus at v = ±4, ±8, ... are visible (dashed lines are a guideline for the eyes). 
The graph in (b) is a plot of the conductance versus gate voltage for B=2T (T=0.3K) showing 
the bilayer Hall quantization sequence. 

Quite recently, the experimental observation of the magnetotransport and the qua ntum Hall 
effect is revealing yet a rich scenario of Landau level spectrum in trilayer graphene iBao et all 



j201l|) : lKumar et alj (|2oTl|) ; lTaychatanapat et alJ (|201l|) ; latang et alj j201l|) l For ABA multi- 
layer graphenes, the Landau level spectra can be again decomposed into a superposition of 
the monolayer and bilayer subsystems as introduced in Sec. |2] In this case, the physical 
properties in magnetic fields, such as Hall conductivity and the magnetic susceptibility can be 
expressed as the summation over components of subsystems iKoshino et al. (2007; 2008)). In 
trilayer graphene, for example, the spectrum is composed of bilayer and monolayer Landau 
levels, resulting in a 12-fold degeneracy at zero energy. The effect of the next-nearest interlayer 
couplings, such as 72 snd 75 (see Fig. |4]a), are often neglected in the simplest approximation, 
but b ecome particularly im portant for the low-energy spectrum near the charge neutrality 
point dKoshino et al.l i201ll) '). For trilayer, the 12-fold degeneracy is lifted by those couplings, 
causing a qualitative change in the quantum Hall plateau structure. The Landau spectrum of 
ABC multilayers is quite di fferent from ABA's . where the pair of low-energy flat-bands gives 
the Landau level sequence jGuinea et all j2006[) ; lKoshino et alj ( 12009.) ): 



Es„=s ^'^^% ^w(w-l)--- (w-N-H), (s = ±, w = 0,1,2,..), (15) 

where N is the number of layers. Including valley and spin degree of freedom, 4N-fold de- 
generate Landau levels appear at zero energy. 

5. Conclusions 

The gate tuneable band structure of FLGs is an unprecedented physical property which paves 
the way to conceptually novel physical phenomena. For instance, an asymmetry induced 
by a perpendicular electric field applied onto bilay er graphene, not only opens a gap bu t 
it also affects the pseudospin of the charge carriers iMin et alj j2008l) ; ISan-Tose et alj feOO^ . 
This pseudospin characterises the layer degree of freedom, and it constitutes an additional 
quantum number for the charge carriers. Recent theoretical schemes propose the use of the 
pseudospin for new devices in which an on/ off state of the current is attained respectively for 
parallel and antiparallel pseudospin configurations in the bilayer. In these pseudospin-valve 
devices the polarity of the electric field acting on bilayer graphene plays a similar role as the 
magnetic field in spin-valve devices. This is the emerging field of pseudo-spintronics. 
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